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ABSTRACT 

We report the discovery of two young isolated radio pulsars with very high inferred magnetic fields. 
PSR Jl 119-6127 has period P = 0.407 s, and the largest period derivative known among radio pulsars, P = 
4.0 X 10"'^. Under standard assumptions these parameters imply a characteristic spin-down age of only = 1 .6 kyr 
and a surface dipole magnetic field strength of Z? = 4.1 x 10^^ G. We have measured a stationary period-second- 
derivative for this pulsar, resulting in a braking index of n = 2.91 ±0.05. We have also observed a glitch in 
the rotation of the pulsar, with fractional period change AP/P = -4.4 x 10"^. Archival radio imaging data sug- 
gest the presence of a previously uncataloged supernova remnant centered on the pulsar The second pulsar, 
PSR J1814-1744, has P = 3.975 s and P = 7.4 x 10"'l These parameters imply = 85 kyr, and B = 5.5 x 10'^ G, 
the largest of any known radio pulsar. 

Both PSR 111 19-6127 and PSR J18 14-1744 show apparently normal radio emission in a regime of magnetic 
field strength where some models predict that no emission should occur. Also, PSR J1814-1744 has spin pa- 
rameters similar to the anomalous X-ray pulsar (AXP) IE 2259H-586, but shows no discernible X-ray emission. 
If AXPs are isolated, high magnetic field neutron stars ("magnetars"), these results suggest that their unusual 
attributes are unlikely to be merely a consequence of their very high inferred magnetic fields. 

Subject headings: pulsars: individual (PSR Jl 1 19-6127, PSR J1814-1744) 



1. INTRODUCTION 

The pulsar in the Crab nebula (PSR B0531H-21), with pe- 
riod f = 33 ms, was born in a type II supernova observed in 
1054 AD, supporting the view that at least some core-collapse 
supernovae (SNe) form pulsars. Based largely on studies of 
the Crab and a few other young objects, a picture has emerged 
where pulsars are born spinning rapidly (with initial period 
Po « 20 ms in the case of the Crab), and spin down due to 
their large magnetic moments according to i> oc -v". In this 
spin-down law = 1/P is the pulsar rotation frequency, i) is its 
derivative, and n = vv/ii))^ is the "braking index." Integration 
of the spin-down law with constant magnetic moment gives the 
age of the pulsar. 



10^ cm and moment of inertia of 10"*^ gcm^, the surface mag- 
netic field strength is 



B = 3.2 
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Braking indices have been measured for 
sars, namely PSRs B0531H-21, B0540-69, 



only four 
B0833-45, 
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The luminosity generated in the braking of the pulsar rotation, 
E = AtP'Ivv, is emitted in the form of magnetic dipole radia- 
tion and a relativistic particle wind. The vast majority of this 
luminosity may be deposited in the ambient environment, pow- 
ering a plerionic supernova remnant (SNR) such as the Crab 
synchrotron nebula, while a very small portion may be observed 
as pulsed electromagnetic radiation. 

Despite the above, many questions remain regarding the out- 
come of type II SNe and the manifestation of young neutron 
stars. Although Galactic SNe and pulsar formation rates are 
both notoriously difficult to estimate (see, e.g., van den Bergh 
& Tam mann 1991; [Tammann, Loffler, & Schroder 1991 Wolt- 



1993 



jer 1998, and tSIarayan & Ostriker 199C| ; l^orimer et al 
Lyne et al. 1998 ), it is quite plausible that type II SNe occur 



B1509-58, with values forw of 2.51 ±0.01, 2.2±0.1, 1.4±0.2, 



and 2.837 ±0.001 respectively (iLyne. Pritchard, & Smith 1993 



Deeter, Nagase, & Boynton 1999|;|Lvne et al. 19961; K aspi et al 



1994). In other cases, an oblique rotating vacuum dipole model 
is ty pically assumed, for which w = 3 (M anchester & Taylor 
1977), and if Pq <C P, equation (|1|) reduces to t = P/{2P) = Tc, 
the characteristic age of a pulsar. With a neutron star radius of 



significantly more o ften than radio pulsars of the kind already 
kno wn are born (see |van den Bergh & Tammann 199l|; W oltj er 
1998). If this is the case, perhaps some young neutron stars are 
being "missed." A possible example is SNR 3C58, the likely 
outcome of a type II SN observed about 820 yr ago, with no de- 
tectable pulsar Studying the energetics and morphology of the 
remnant, Helfand, Becker, & White (1995) make a compelling 
case for the presence of an unseen pulsar with higher magnetic 
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field than any previously known. Having a short period like the 
Crab at birth, such a pulsar would have spun down rapidly to 
a present long period. Maybe yet other pulsars are born spin- 
ning slowly and never generate the large E required to power an 
easily detectable nebula. In addition, some neutron stars may 
never manifest themselves as radio pulsars at all. It has been 
suggested that there exists a class of isolated rotating neutron 
stars with ultra-strong magnetic fi elds, the so-called "magne- 
tars" ( puncan & Thompson 1992 ). The observational proper- 
ties of radio pulsars and magnetar candidates are very different. 
Radio pulsars rarely exhibit X-ray pulsations, and when they 
do, their X-ray power is small compared to their £. By contrast, 
magnetars emit pulsed X-rays with luminosities fa r in excess of 
their spin- down power (|Vasisht & Gotthelf 1997|; K ouveliotou 
et al. 1999) but remain undetected at radio wavelengths. The 
dichotomy is thoug ht to result from the much l arger magnetic 
fields in m agnetars ([Thompson & Duncan 1993t H eyl & Hern- 
quist 1997). 

In this paper we report the discovery of two isolated radio 
pulsars with some properties that are unusual and interesting in 
the context of the above questions. 

2. OBSERVATIONS AND RESULTS 

The radio pulsars J1119-6127 and J1814-1744 were dis- 
covered on 1997 August 24 and 23, respectively, in a survey 
of the Galactic plane using the 64-m Parkes r adio telescope 
in A ustralia. This survey (^ne et al. 2000a|; C amilo et al. 
2000) makes use of the fast rate of sky coverage afforded by 
a multibeam receiver to increase greatly the integration time, 
and consequently the sensitivity, relative to previous surveys. 
The Parkes survey uses 13 beams at a central sky frequency of 
1374 MHz with an equivalent system noise of = 35 Jy at high 
Galactic latitude. For each beam, the sum of two orthogonal lin- 
ear polarization channels, each 288 MHz wide, is recorded for 
35 minutes per grid position, providing sensitivity to all pulsars 
with flux densities in excess of ^ 0.15 mJy for P ^ 0.1 s. 

Follow-on regular timing observations have been carried out 
at Parkes since 1998 February for newly discovered pulsars 
with declination south of -35°, while the remainder are ob- 
served with the 76-m telescope at Jodrell Bank Observatory, 
England. The system used for timing observations at Parkes 
is identical to that used in the survey, although we record sig- 
nals from the central beam only: the down-converted radio- 
frequency noise is passed through a 2 x 96 x 3-MHz filter bank 
spectrometer, after which the signals are square-law detected, 
orthogonal polarizations are summed, and the 96 resulting volt- 
ages are high-pass filtered before being 1-bit digitized every 
250 /is and written to magnetic tape for subsequent analysis. 
We also record the start time of each observation, synchro- 
nized with the observatory time standard and traceable to UTC. 
PSR Jl 1 19-6127, in whose direction 5o = 40 Jy, was observed 
in this manner on 63 days over a period of two years, for 
approximately 10 minutes each day. At Jodrell Bank, with 
Sq = 50 Jy in the direction of PSR J18 14-1744, the observing 
setup used a 2 x 32 x 3-MHz filter bank until 1999 July, and a 
2 X 64 X 1-MHz filter bank since then, to observe a band cen- 
tered in the range 1376 to 1396 MHz, depending on the radio- 
frequency interference environment. Signals from individual 
frequency channels are delayed by an amount proportional to 
the dispersion measure (DM) of the pulsar, to account for dis- 

'http://pulsar.princeton.edu/tempo. 
'"http ://w w w. atnf .csiro . au/computing/ software/miriad. 



persion caused by propagation through the interstellar medium, 
and are folded synchronously with the predicted rotation pe- 
riod, generating one pulse profile for each sub-integration last- 
ing 3 minutes. We have observed PSR J1814-1744in this man- 
ner on 37 days over a two-year interval, for 18 minutes each 
day. 

We have analyzed the timing data in standard fashion. 
Briefly, topocentric pulse times-of-arrival (TOAs) were mea- 
sured by cross-correlating daily-averaged pulse profiles with 
a high signal-to-noise-ratio template (see Fig. p. Celestial 
coordinates and spin parameters were then determined using 
the TEMPO software pack age^ and the JPL DE200 planetary 
ephemeris ( [Standish 199C ). TEMPO first converts the TOAs to 
the barycenter of the solar system, and refines the initial esti- 
mated parameters in a fitting procedure that minimizes timing 
residuals (difference between observed and predicted TOAs) 
with respect to the model parameters. 





Pulse phase (sec) 

Fig. 1.— Integrated pulse profiles for PSRs Jl 119-6127 (top) and 
J1814-1744 (bottom) at a frequency of 1374 MHz. 

Underlying the timing model is the assumption that the rota- 
tional phase of the neutron star is described by 



(3) 



1 , 1 

2 6 

where T denotes pulsar proper time. In equation (|3]), the inter- 
pretation of ly and its derivatives as representing only the sta- 
tionary spin parameters of a rotating magnetic dipole does not 
hold strictly for pulsars displaying rotational irregularities, and 
parameter estimation in such circumstances must be performed 
with extra care. 

PSR Jl 119-6127 was observed on 1998 October 30 and 31 
with the Australia Telescope Compact A rray (ATCA), wi th the 
interferometer in its "6D" configuration (Lazendic 1999). The 



observations were done in pulsar gating mode simultaneously 
at center frequencies of 1384 and 2496 MHz, with 128 MHz of 
bandwidth in each of two linear polarizations at each frequency. 
The radio sources 1934-638 and 1036-697 were used as flux 
density and phase calibrators, respectively. The data were pro- 
cessed using the MIRIAD package'", during which on- and off- 
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pulse maps were generated. The data set at 1384 MHz, at which 
frequency the pulsar is brightest, was used to obtain the posi- 
tion of the pulsar, and the 2496 MHz observation was used to 
determine the flux density, both listed in Table |l} 

Because PSR J 111 9-6 127 has the largest period derivative 
(P) of any kno wn radio pulsar, and P co rrelates highly with 
"timing noise" ( Arzoumanian et al. 1994) which can bias the 



the residuals. Arzoumanian et al. (1994) use the parameter 

1 



A(f) = iog —\iy\t 



(4) 



celestial coordinates determined with timing data, we use the 
position obtained from interferometric observations in the tim- 
ing solution. We then fit for v, i>, and v. Figures ^ and 
b indicate a small glitch in rotation occurred on about MJD 
51398; glitch parameters are given in Table [l|. Figure ^ also 
suggests that another glitch of similar magnitude may have oc- 
curred sometime during MJD 50850-50940, but we cannot be 
sure. The rotational parameters best describing the behavior of 
the pulsar are listed in Table |l[ with the corresponding timing 
residuals shown in Figure 0c. 
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Fig. 2. — (a) Timing residuals for PSR Jl 1 19-6127 using data collected 
over 1998 May-1999 August after fitting for u, u, and u. (b) Residuals for en- 
tire data set using model obtained in (a), showing a glitch at about MJD 51398. 
(c) Residuals using data collected over 1998 May-20()0 February after fitting 
for u, i), V, and glitch parameters, (d) Post-fit residuals for PSR J1814-1744. 



In Figure ^ the residuals following the glitch appear cubic 
in shape, with amplitude much reduced by comparison with the 
parabolic residuals in Figure ^7. This suggests that the glitch 
parameters in Table [l| do not completely describe the behavior 
of the pulsar following the event. Further data are required to 
determine whether the post-glitch spin parameters are relaxing 
with exponential decay time-scales of order several months, as 
seen in the Vela pulsar, or whether the change in v at the glitch 
(or at least much of it) is permanent, as observed in the Crab 
pulsar (see Lyne, Shemar, & Graham-Smith 2000b| ). 

The uncertainty in the braking index determined from the 
spin parameters, n = 2.91 ±0.01, reflects only the random phase 
noise resulting from uncertainty in the TOAs. We now consider 
the effect upon this measurement of the possible presence of 
timing noise and the known occurrence of glitch(es). For most 
pulsars the stationary value of v in equation (|^) is too small to 
be measured, and in timing fits where only v and v are deter- 
mined, excess noise manifests itself as a quasi-cubic trend in 



to estimate the cumulative phase contribution over time t due 
to timing noise and find that, for t= 10^ s (arbitrary, but similar 
to the time span of their observations), most pulsars, despite a 
large scatter in the data, follow the relationship 

A8 = 6.6 + 0.61ogP. (5) 

Using equations (§)-(§) for PSR Jl 1 19-6127, we estimate that 
the measured value of v may be contaminated by as much 
as 8 X 10~^'*s"-', or four times the formal uncertainty given in 
Table |l]. In fact we have not measured any timing noise for 
this pulsar, as indicated by the apparently "white" residuals 
(Fig. ^) and by the upper limit on, rather than measurement 
of, V obtained in an additional fit to the data represented in Fig- 
ure |fl (see Table [l]). However this is not too surprising, given 
that our inter-glitch data span only 1 .2 yr; with a longer time 
span between glitches, timing noise, if it is present, may be 
measurable. In summary, we believe an accurate measurement 
of the braking index between glitches, reflecting the steady 
spin-down physics of the neutron star, is n = 2.91 ±0.01 ±0.04. 

If the change in i> at the glitch, Ai' (Table |l]), is perma- 
nent, it contributes a component to v beyond that measured 
between glitches. This contribution is approximately At>/Ar, 
where AT is the time interval between glitches. By assum- 
ing that all of the measured Ai> is permanent, and that a first 
glitch did occur at MJD ^ 50900 (see Fig. ^), we obtain 
Aj>' = -2 X 10"^^ s"-', about 10 times the formal uncertainty in 
V and implying a correction to n of -0.1. Future measurements 
will settle this question, but for the purposes of calculating pul- 
sar age, the correct value of braking index may be as low as 
« 2.8. A similar effect is seen in the Crab pulsar, where perma- 
nent changes in spin-down rate occurring at glitches contribute 
a correct ion of -0.0 5 to the value of n = 2.5 measured between 
glitches ( [Nice 1993[ ). 

The quasi-cubic trend in the residuals of PSR J1814-1744 
shown in Figure ^ suggests the presence of timing noise in 
this pulsar. In these circumstances, we determined its position 
by "whitening" the residuals with a fit of the data to a model 
incorporating celestial coordinates, v, i), and v. We took the 
resulting coordinates, with uncertainties, as our best unbiased 
estimate of position, and kept them fixed in subsequent fits to 
the timing data. We then performed one fit for v and i), result- 
ing in the best values for these stationary parameters averaged 
over the data span, with residuals displayed in Figure ^ and 
showing some red-noise. Finally we performed one extra fit 
with the additional free parameter v, which is not stationary. 
Values obtained for all these parameters are listed in Table nl 
We note that the value of v, an estimate of the amount of timing 
noise present in PSR J1814-1744, is approximately at the level 
expected from equations (Q)-(^). Note that the uncertainty in 
declination is particularly large because this pulsar is located at 
low ecliptic latitude. 

The DMs quoted in Table |l] were obtained by folding raw 
data at the known pulse period for each of four frequency 
sub-bands, created by the addition of data from 24 adja- 
cent frequency channels, and fitting a time-delay between the 
sub-bands. This was done for timing data in the case of 
PSR J 1 1 1 9-6 1 27 and discovery data for PSR J 1 8 14- 1 744. Dis- 
persion measures, togethe r with a model for the G alactic elec- 
tron density distribution (Taylor & Cordes 1993), are used to 



4 



Table 1 

Measured and derived parameters for pulsars Jl 1 19-6127 and J1814-1744 



Right ascension, a (J2000) 

Dechnation, S (J2000) 

Rotation frequency, z-* (s~' ) 

Frequency derivative, (s~^) 

Second frequency derivative, (s~' ) 
Third frequency derivative, (s""*) . 

Epoch (MJD) 

Frequency step at glitch, Ai/ (s~^) . 
Change in at glitch, Ai> (s~^) .... 

Epoch of glitch (MJD) 

R.M.S. residual (ms) (white/red) . . . 
Dispersion measure, DM {cm~^ pc) . 
Flux density at 1374 MHz, 5 (mJy) . 
Flux density at 2496 MHz (mJy) . . . 

Spin period, P (s) 

Period derivative, P 

Surface magnetic field, B (Gauss) . . 

Characteristic age, (kyr) 

Spin-down luminosity, E (ergs"') . . 

Braking index, n 

Distance, d (kpc) 

Radio luminosity, Sd~ (mJykpc^) . . 

Galactic longitude, / (deg) 

Galactic latitude, b (deg) 



PSRJl 119-6127 


PSR J1814-1744 


11 1 n 1 /I '3a/i\l| 
illy 14. iU(Z) 


10 1/1 /1TA/1/1A\ 

18 14 4Z. 94(10) 


/^l 01 AQ ^("01'' 


— 1 / ZJ^ ly) 


O /iClI*^A1 I 1 ^ 
Z.4j J lOUl 1 JU(i ) 


U.Zj 1 J ly /41 J^J ) 


-2.4207996(8) x 10"" 


-4.7002(4) X 10"' 


6.94(2) X 10"^^ 


3.0(7) X 10"^" '' 


< 10-" 




51173.0 


51200.0 


1.08(10) X 10"* 




-9.5(13) X ur" 




51398(4) 




0.7/ • • ■ 


6.3/11.2 


707(2) 


834(20) 


0.9(1) 


0.8(1) 


0.44(5)" 




0.40763747736(2) 


3.975831061(5) 


4.022602(2) x 10"'^ 


7.4297(6) X 10"" 


4.1 X 10" 


5.5 X 10" 


1.6 


85 


2.3 X 10^" 


4.7 X 10^^ 


2.91(1) 




2.4-8 


10(2) 


~ 25 


~ 80 


292.15 


13.02 


-0.54 


-0.21 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. Figures in parentheses represent 1 a unceitainties in 
least-significant digits quoted. 

^Obtained from interferometric observations (see Celestial coordinates obtained from a fit to data collected over MJD 5094()-5 1392 are a = ll''19°'14!24(5), i5 = -61°27'49"8(5). 
^This parameter is not stationary (see § ^ for details). 



estimate distances to pulsars. For PSR Jl 1 19-6127 the implied 
distance is li > 30 kpc. However, the model does not account 
for most individual Hll regions, and, because the pulsar lies in 
the direction of the Carina spiral arm, with the likelihood of 
nearby ionizing population I stars, this distance is assuredly a 
gross overestimate. We believe it is likely that the pulsar is lo- 
cated between the two line-of-sight crossings of the Carina arm, 
between d = 2A and 8 kpc. The distance estimate quoted in Ta- 
ble IT] for PSR J18 14-1744 is that obtained from the model of 
Taylor &Cordes (1993). 

The flux densities at 1374 MHz listed in Table |l] were de- 
termined by converting the average observed signal strength of 
the pulsars to a scale calibrated using published flux densities at 
1400 MHz for a group of high-DM pulsars, taking into account 
the variation in sky background temperature. 

3. DISCUSSION 

3.1. PSR Jl 119-6127 

PSR Jl 119-6127 has the largest period derivative known 
among radio pulsars. Partly for this reason it was rela- 
tively straightforward to measure a stationary i/ with a phase- 
connected timing solution (i.e., through absolute pulse number- 
ing), only the third pulsar for which this has been possible. The 
resulting value of braking index is « = 2.91 ±0.05, including 
possible contamination by timing noise (see § and is in good 
agreement with that predicted by a model treating the pulsar as 
an oblique rota tor with a current-starved outer magnetosphere 
(Melatos 1997). For the four other pulsars for which it has been 
measured, n ranges between 1.4 and 2.8 (see § |l]). That ob- 



served braking indic es are smaller than 3 can be explained in a 
variety of ways (see Melatos 1997 for a review), including a ki- 
netic energy-dominated flow at the light cylinder, or an increase 
in the mag netic moment of the star over time (B landford & Ro- 
mani 1988). None of these scenarios are consistent with all the 



observations ( Arons 1992). Measurement of iy would constrain 
these possibilities further. At present the upper limit in Table [l| 
i s 30 times the value expect ed from a simple spin-down law 
(Blandford & Romani 1988). Whether this can be measured, 
and how much the measurement of n can be improved with fur- 
ther observations, will depend on the level of timing noise and 
glitch activity displayed by the pulsar. 

Assuming that Pq <C P, but using the measured values of P, 
P, and n (Table [^) in equation (|T|), the age of PSR Jl 1 19-6127 
is 1.7 ±0.1 kyr, including possible biases due to timing noise 
and glitches (see § Of course, if the pulsar were born 
spinning slower, it would be younger. For Pq = 0.2 s, half the 
present period, the age is 1.2 kyr. In any case it is clear that 
PSR Jl 119-6127 is among the very youngest neutron stars 
known. 

Three other pulsars with characteristic ages under 2 kyr are 
known: the Crab pulsar (t^ = 1.3 kyr), PSR B1509-58 in 
G320.4-1.2 (tc = 1.6 kyr), and PSR B0540-69 in the Large 
Magellanic Cloud (t^ = 1.7 kyr). All three are associated 
with SNRs. We have searched for evidence of an SNR nea r 
PSR Jl 119-6127. Although none is cataloged ( preen 1996| ), 
data from the Molonglo Observatory Sy nthesis Telescope ob- 



tained at a radio frequency of 843 MHz ( Green et al. 1999| ) re- 
veal a faint ring of radius 7' centered on the pulsar. This could 
be the expanding blast wave of the SNR. Its size would imply an 
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expansion velocity of lO'* km s~', for a distance of 8 kpc and 
age of 1 .6 kyr, reasonable for the blast- wave interpretation if the 
surrounding medium is of low density and relatively uniform. 
Additional ATCA data show that the shell has a non-thermal 
radio spectrum (Crawford et al., in preparation). This possible 
SNR is also X-ray-bright, with its spectrum described by either 
a power-law or thermal model, and additional observations are 
required to furth er constrain its properties (P ivovaroff, Kaspi, 
& Camilo 2000a). Although the supernova that gave birth to 
this pulsar occurred in an era in which celestial events were 
recorded by some civilizations, this explosion may have been 
too far south and/or too distant or too obscured to have been 
detected by these observers. 

The glitch observed in PSR Jl 119-6127 is small compared 
to most glitches in most pulsars, with Aj//i/ = (4.4±0.4) x 10"^ 
(Table |]), but it is of similar fractional size as three of the five 
glitches observed in the Crab pulsar over 23 years (Lyne et 
al. 2000b). It remains to be seen whether at least some of the 
change measured in v is permanent, as seen in the Crab glitches. 
Unless we were unreasonably lucky, PSR J 11 V^-dVLl glitches 
more often than the Crab pulsar, but it is curious that its glitches 
share some characteristics with those of the Crab: while its pe- 
riod and magnetic field are approximately 10 times larger than 
the Crab's, its age, and perhaps therefore its internal tempera- 
ture, are similar 

Finally, we compare the PSR Jl 1 19-6127 system with some 
young pulsar/SNR systems. For ages <^ 2000 yr, the radio lu- 
minosity Lr of a synchrotron nebula ("plerion") with a cen- 
tral pulsar is a measure of the energy output of the pulsar 
over its lifetime, due to the relatively long lifetime of the ra- 
diating electrons. The plerion X-ray luminosity Lx, on the 
other hand, reflects the current E of the pulsar. For the Crab 
and PSR B0540-69, Lx 0.05£, while for PSR B1509-58, 
Lx ^ 0.01/i (see Helfand et al. 1995, and references therein). 
For SNR 3C58, Helfand et al. find that all available data can be 
reconciled with a (candidate) pulsar having P ^ 0.2 s, P 4 x 
10"'^ (parameters similar to those of PSR Jl 119-6127 — see 
Table |l]), and with Lx - 5 x \<ir^E. For PSR Jl 1 19-6127, with 
a current E 200 times smaller than the Crab's, the limit on pleri- 
onic X-ray emission is Lx ^ lO'^L (Pivovaroff et al. 2000a). If 
PSR Jl 1 19-6127 were born with a small period, it would have 
had a much larger E within the past ~ 1700yr, possibly larger 
than the Crab's initially. That energetic past might be reflected 
in plerionic radio emission near the pulsar, depending on the lo- 
cal environment. A measurement of Lx and Lr may in principle 
provide information about whether PSR Jl 1 19-6127 was born 
with a rapid spin rate, as commonly assumed for most pulsars, 
or whether it was born a slow rotator 

3.2. PSRJ1814-1744 

Figure || is a plot of P versus P for the radio pulsar popula- 
tion. PSRs Jl 1 19-6127 and J1814-1744 are indicated, and we 
infer B = 4.1 x lO'-' and 5.5 x lO'-' G respectively, using equa- 
tion (^. These are the highest magnetic field strengths yet ob- 
served among radio pulsars. The pulsars with the next largest 
values of B are PSRs J1726-3530 (P = l.ls;B = 3.7 x lO'^G) 
and J1632-4818 (P = 0.8s;B = 2.3 x lO'^G), also discovered 
in the multibeam survey". Prior to this survey the largest value 
was B = 2.1 X 10'^ G for the 2.4s PSR B0154-H61 (A rzouma- 
nian et al. 1994). 

Also shown in Figure |^ are the sources usually identified as 



magnetars, namely the five anomalous X-ray pulsars (AXPs) 
and two soft gamma repeaters (SGRs) for which P and P have 
been measured. AXPs are characterized by X-ray periods in the 
rang e 5-12 s and extremely rapid spin down (G otthelf & Vasisht 
1998), while the SGRs exhibit occasional enormous bursts of 
7-radiation and AXP-like X-ray pulsations during quiescence. 

Most mo dels of the radio emission physics (M anchester & 
Taylor 1977) depend on pair-production cascades above the 
magnetic poles and hence on the magnitude of the magnetic 
field. However, at field strengths near or above the so-called 
"quantum critical field," 



eh 



:4.4x 10'^ G 



(6) 



the field at which the cyclotron energy is equal to the electron 
rest-mass energy, processes such as photon splitting may inhibit 
pair-producing cascades. It has therefore been argued (B aring 
& Harding 1998) that a radio-loud/radio-quiet boundary can be 
drawn on the P-P diagram, with radio pulsars on one side, and 
AXPs and SGRs on the other (see dotted line in Fig. S). 
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Fig. 3. — Plot of P versus P for radio pulsars (dots), anomalous X-ray pul- 
sars (AXPs), and soft gamma-ray repeaters (SGRs). PSRs Jl 119-6127 and 
J1814-1744 are identified by large filled circles, and sources plausibly associ- 
ated with supernova remnants (SNRs) are noted. Lines of constant character- 
istic age and surface magnetic field strength are drawn. The dotted line shown 
between the lines for B = 10'^ and 10'** Tt indicate s a hypothesized approxi- 
mate theoretical boundary ( Baring & Harding 1998 ) separating radio-loud and 
radio-quiet neutron stars diie to effects relatmg to magnetic fields close to the 
critical field (see discussion following equation [™). 



The existence of PSRs Jl 119-6127, J1726-3530, and 
J1814-1744 demonstrates that radio emission can be produced 
in neutron stars with B^B^. The radio luminosities of these ob- 
jects (Table 1) are typical for observed radio pulsars. Thus, pho- 
ton splitting does not appear to inhibit radio emission at these 
magnetic fields, in agreement with Usov & Mekose (1995) who 
argue that this process is inhibited by polarization selection 
rules. Also, there are both astrophysical and instrumental selec- 
tion effects which bias searches against the detection of long- 
period (P > 5 s) radio pulsars such as Jl 8 14-1744: evidence 
suggests their beams are narrower (e.g.. Young, Manchester, & 



' ' See http://www.atnf.csiro.au/~pulsar/psr/pmsurv/pmwww/pmpsrs.db. 
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Johnston 1999), so the chances of one intersecting our Hne-of- 
sight are smaller, and instrumental high-pass filtering intended 
to remove baseline variations reduces the sensitivity of searches 
for such pulsars. Pulsars such as J1814-1744 could therefore 
be more prevalent than present numbers suggest. 

Especially noteworthy is the proximity of PSR J1814-1744 
to the cluster of AXPs and SGRs at the upper right corner of 
Figure ^ In particular, this pulsar ha s a very similar P to that of 



the well-known AX P IE 2259H-586 ( Fahlman & Gregory 1981 



Baykal et al. 1998), which has a period of 7 s. The disparity in 



their emission properties is therefore surprising. 

The absence of X-ray emission from the direction of 
PSR J1814-1744, inferred from archival ASCA and ROSAT 
observations, impli es that it must be significantly less lum inous 
than IE 2259H-586 ( [Pivovaroff, Kaspi, & Camilo 2000b[ ). 

The radio e mission upper limit for IE 2259-H586 (C oe, Jones, 
& Lehto 1994; |Lorimer, Lyne, & Camilo 1998 ) implies an up- 
per limit on the radio luminosity at 1400 MHz of 0.8 mJy kpc^, 
10"^ that of PSR J1 814-1744, assuming a distance of 4kpc 
(Rho & Petre 1997). This limit is comparable to the lowest 



valu es observed for the radio pulsar population ( Tauris et al. 
1994). That the radio pulse may be unobservable because of 
beaming cannot of course be ruled out. 

The radio-loud/radio-quiet boundary line displayed in Fig- 
ure is more illustrative than quantitative ( Baring & Harding 
1998). However, the apparently normal radio emission from 
PSRs Jl 119-6127, J1726-3530, and J1814-1744, and the ab- 
sence of radio emission from IE 2259H-586, located very close 
to PSR J1814-1744 on a P-P diagram (Fig. H), suggests that it 
may be difficult to delineate any such boundary. 

The two sources are also similar in their levels of rotational 
stabiHty, at least on time scales of ^ 2yr: PSR J1814-1744 
displays timing noise in the amount expected for a radio pul- 
sar with its P (see § as is the upper limit on timing noise 
for IE 2259-H58 6 over a 2-3 yr span (K aspi, Chakrabarty, & 
Steinberger 1999). However, longer term incoherent timing of 
IE 2259H-586 has revealed significant deviations from a sim- 
ple spin-down model. These have been interpreted as being 
evidence for radiative precession of the neutro n star, due to it s 
physical distortion by the strong magnetic field (Melatos 1999). 
Alternatively, Heyl & Hernquist (1999) suggest the deviations 
are due to extremely large glitches. In either model, similar be- 



havior might be expected of PSR J1814-1744; continued radio 
timing will be sensitive to it. 

The similar spin parameters for these two stars and, in turn, 
many common features between IE 2259H-586 and some other 
AXPs and SGRs, suggest that very high inferred magnetic field 
strengths cannot be the sole factor governing whether or not an 
isolated neutron star is a magnetar or a radio pulsar Other pos- 
sible facto rs include heavy-element a t mospheric composition 
and youth ( Thompson & Duncan 1 993 ; [Heyl & Hernquist 1997 



see also Pivovaroff et al. 2000b). The age of PSR J1814-1744, 
if P[) ^ P and « = 3, is 85 kyr It is unlikely that any associated 
supernova remnant would stil l be observab le and indeed there 
is none known in the vicinity (Green 1996). 

We also note that the recently proposed accret ion model 
for AXPs (Chatterjee, Hernquist, & Narayan 2000), in which 
they are accreting from a fall-back disk formed from mate- 
rial remaining after the supernova explosion, is challenged by 
PSR J1814-1744. In this model, the neutron star should not 
be a radio pulsar, but rather an AXP progenitor in a "dim pro- 
peller phase," its rotational frequency being still too high for 
the accreting material to overcome the centrifugal barrier Of 
course, it is always possible that in this one case no fall-back 
disk formed. 

Proof that AXPs or SGRs are isolated high-magnetic- 
field neutron stars would come from either the discovery of 
magnetar-like emission from a radio pulsar, or radio pulsations 
from a putative magnetar While such radio emission was not 
expected due to theoretical considerations, because of the high 
inferred magnetic fields, the discovery of PSR J1814-1744 
shows that this emission does occur at magnetic field values 
characteristic of at least some magnetars, opening the possibil- 
ity that magnetars also emit observable radio waves. 
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